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Abstract. We have studied the problem of phase stability in NiPt alloy system. 
We have used the augmented space recursion based on the TB-LMTO as the method 
for studying the electronic structure of the alloys. In particular, we have used the 
relativistic generalization of our earlier technique. We note that, in order to predict 
the proper ground state structures and energetics, in addition to relativistic effects, we 
have to take into account charge transfer effects with precision. 

PACS numbers: 71.20,71.20c 
1. Introduction. 

There has been growing interest in the study of alloy phase ordering and segregation 
using first principles techniques. In order to study these phenomena one needs a 
derivation of the configurational energy for the alloy system. Different models have 
been proposed in which the configurational energies are expressed in terms of effective 
multi-site interactions, in particular effective pair interactions [1. The analysis of 
alloy ordering tendencies and phase stability reduces to the accurate and reliable 
determination of effective pair interactions . There are two different approaches of 
obtaining the effective pair interactions. One approach is to start with electronic 
structure calculations of the total energy of ordered super-structures of the alloy and to 
invert these total energies to obtain the effective pair interactions. This is the Connolly- 
Williams method [2]. The other approach is to start from the completely disordered 
phase, set up a perturbation in the form of concentration fluctuations associated with 
an ordered phase and study whether the alloy can sustain such a perturbation. This 
includes approaches like the generalized perturbation method (GPM) 3J, the embedded 
cluster method (ECM)[4J. Most of the works on the electronic structure of the disordered 
alloys have been based so far on the coherent potential approximation (CPA). The CPA 
being a single-site approximation cannot take into account the effect at a site of its 
immediate environment. In an attempt to go beyond the single site approximation, 
de Fontaine and his group followed a different approach of direct configurational 
averaging (DCA) jH], without resorting to any kind of single-site approximation. The 
effective pair and multi-site interactions were calculated directly in real space for given 
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configurations and the averaging was done in a brute force way by summing over different 
configurations. Invariably, the number of configurations was finite and convergence of 
the results with increasing number of configurations is yet to be available. 

Saha et al |B| have introduced the augmented space recursion (ASR) based on the 
augmented space formalism (ASF) first suggested by Mookerjee [7j coupled with the 
recursion method of Haydock et al jSj- Within ASF the configuration averaging is 
carried out without having to resort to any single-site approximation. The recursion 
method allows us to take into account the effect of the environment of a given site . 
Moreover, the convergence of various physical quantities calculated through recursion 
with the number of recursion steps and subsequent termination has been studied in great 
detail [Oj. Among advantages of the ASR in going beyond the single-site approximation 
is the possibility of inclusion of local lattice distortions which is important in the case 
of alloys with large size mismatch between components as in the case of NiPt. In an 
earlier paper Saha and Mookerjee jTU] had discussed the effect of local lattice distortion 
on the electronic structures of CuPd and CuBe alloys using the ASR. This allows the 
structure matrices to randomly take values Sif,, S^E'i Sjtv or ^El'j depending on the 
occupation of the sites R and R'. The augmented space recursion coupled with orbital 
peeling technique [TT] to evaluate small energy differences associated with band structure 
energies has been successfully used in past to describe the phase formation in alloys [T2] . 

In the present communication we focus on the application of this method for phase 
stability study in NiPt alloys. This system of alloys is of importance because of the 
possible need for relativistic corrections due to the heavy mass of Pt as well as effects 
due to charge transfer and size mismatch between Ni and Pt. This therefore forms a 
perfect candidate for testing the applicability and limitations of our formalism, bringing 
in the relative importance of various effects for the accurate description of the system. 
The previous studies of ordered and substitutionally disordered NiPt alloy systems have 
shown the importance of inclusion of relativistic effects. Treglia and Ducastelle J3] 
had shown that late transition metal alloys should exhibit phase separating tendencies 
but they argue that the exceptional ordering behavior of NiPt is due to the relativistic 
corrections. In a first principle study, Pinski et al [13] found that the disordered fee 
Nii_ x Pt x alloy at x — 0.5, calculated by means of the single site KKR-CPA, becomes 
unstable at low temperatures, to a perturbation by a (100) ordering wave and concluded 
that the corresponding long range ordered state (LRO) i.e. the Ll structure should be 
the predicted ground state for which the large size mismatch between Ni and Pt plays the 
main role and the effect of relativity can be neglected. However, Lu et al [15J pointed out 
that a local ordering tendency determined by perturbation analysis, doesn't necessarily 
predict the correct LRO ground state if the size mismatch of the two elements is large, 
as is the case for Ni and Pt and concluded that relativity is the sole reason for long range 
order in NiPt. The work of Singh et al jTHj demonstrated that the relativistic effects do 
stabilize the ordered structures over the disordered solid solution. Recently Ruban et al 
[T7] have studied the problem of phase stability in NiPt alloy system based on ordered 
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calculations with the inclusion of Madelung energy with multipole corrections. In this 
paper, we examine the relativistic treatment of the Hamiltonian and charge transfer 
and lattice relaxation effects on the electronic structure and phase stability of face- 
centered cubic NiPt system at 25%, 50% and 75% of concentration of Pt. As mentioned 
already, the augmented space recursion (ASR) technique, which we use here, is capable 
of taking into account environmental effects, effects of short range order and local lattice 
relaxation effects due to size mismatches. To circumvent the problem of calculation of 
Madelung energy contribution for disordered system, we have used the appropriate 
effective atomic sphere radii for each of the constituents so that the spheres are neutral 
on the average and this has been done with precision at each concentration ^Hj- We 
have shown that without inclusion of relativistic effects the formation energy comes 
out to be positive which contradicts experimental results. With the scalar relativistic 
corrections, involving mass-velocity and Darwin terms, the formation energy comes out 
negative indicating that the relativistic effects play an important role in NiPt alloys in 
agreement with earlier studies. We find that the charge transfer effects have also an 
important role to play in deciding on the correct ground state structure, particularly 
when the concentration of Pt is high. Our study on transition temperatures based on a 
mean field theory could reproduce the qualitative experimental trends. 

2. Formalism 

2.1. The Effective Pair Interactions 

We start from a completely disordered alloy. Each site R has an occupation variable 
n R associated with it. For a homogeneous perfect disorder (%) = x, where x is 
the concentration of one of the components of the alloy. In this homogeneously 
disordered system we now introduce fluctuations in the occupation variable at each 
site : 5xr = Ur — x. Expanding the total energy in this configuration about the energy 
of the perfectly disordered state we get : 

N N 

E{x) = + J2 E R 5x R + E E RR> 5x R 5x R' + ■■■ (!) 

R=l RR'=1 

The coefficients £ (0) , E$ . . . are the effective renormalized cluster interactions. E^ 
is the energy of the averaged disordered medium. The renormalized pair interactions 
E RRI express the correlation between two sites and are the most dominant quantities 
for the analysis of phase stability. We will retain terms up to pair interactions in the 
configuration energy expansion. Higher order interactions may be included for a more 
accurate and complete description. For the phase stability study it is the pair interaction 
which plays the dominant role. 

The total energy of a solid may be separated into two terms : a one-electron band 
contribution Ebs and the electrostatic contribution Ees The renormalized cluster 
interactions defined in (|TJ) should, in principle, include both E BS and E ES contributions. 
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Since the renormalized cluster interactions involve the difference of cluster energies, it 
is usually assumed that the electrostatic terms cancel out and only the band structure 
contribution is important. Such an assumption which is not rigorously true, has been 
shown to be approximately valid in a number of alloy systems [19|. Considering only 
band structure contribution, the effective pair interactions may be written as : 

4% = -jZ dE {~^ m lo ^det(G IJ )(E)^ (2) 

where, G IJ represents the configurationally averaged Green function corresponding to 
the disordered Hamiltonian whose R and R' sites are occupied by I-th and J-th type of 
atom, and 



f +1 if I=J 
\ -1 if M J 

The behavior of this function is quite complicated and hence the integration by standard 
routines (e.g. Simpson's rule or Chebyshev polynomials) is difficult, involving many 
iterations before convergence is achieved. Furthermore the integrand is multi-valued, 
being simply the phase of J2ij det (G /J ) The way out for this was suggested 

by Burke |lljwhich relies on the repeated application of the partition theorem on the 
Hamiltonian H /J . The final result is given simply in terms of the zeroes and poles of 
the Green function in the region E < Ep 



Imax 

E RR' = 2 J2&J E 
IJ k=0 

rkJJ j r>k,IJ 



pk.u 



E Z V J - E P V J + ( P kJJ - z k > IJ ) E t 



(3) 



where Z, ' and P- ' are the zeros and poles of the peeled Green's function G{ J of 
disordered Hamiltonian with occupancy at sites R and R' by I and J of which first (k-1) 
rows and columns has been deleted. p h ' IJ and z k,IJ are the number of poles and zeroes 
in the energy region below Ep. The factor 2 accounts for the spin degeneracy. 



2.2. The Augmented Space Recursion 

As discussed in the previous section, the calculation of the effective pair interaction in 
our formalism reduces to the determination of the peeled configuration averaged green 
functions (G^ J ). We shall employ the augmented space recursion coupled with the 
linearized tight-binding muffin tin orbital method (TB-LMTO) introduced by Andersen 
and Jepsen j20j for a first principles determination of these configuration averaged 
quantities. We shall take the most localized, sparse tight binding first order Hamiltonian 
derived systematically from the LMTO theory within the atomic sphere approximation 
(ASA) and generalized to random alloys. The augmented space recursion method has 
been described at great length in earlier communications ([Hl-IH], ^2] , EI] , [22] ) • We refer 
the readers to these references for the details. We shall give here the final form of the 
effective Hamiltonian used for recursion in augmented space for the calculation of the 
peeled Green functions : 
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Hl J = E Ci/^n + £ o^a w + E E + 5CM R ") 4,a*' + ... 

t=k 1=1 R"^R,R' i=l 

+ E EEA^fi?''(A^ + 5AfM-'') aW + - 
+ E E E K?f + ^ 2 M«") A^¥ aU* + . . . 

R"^R L L'=k 

+ E EE ( A ii'V + 6A i M J LL' A R>j>> a R» a R> + ■■■ 

R"jtR' L L> 

+ £ E EE (Aff + iAfM")*' (A^+<M'')... 

R"^R,R' R"'^R,R' L L' 

■ ■ ■ {a ] R „a R n, + a Rlll a R „) (4) 

here, L is a composite index (Im) 

For a binary distribution is given by: 

M R = x b R] b m + (1 - x) b Rl b Rl + y/x(l-x) (b R] b Rl + 6^6^) (5) 

For non-isochoric alloys , the difference in atomic radii of the constituents lead to change 
in the electronic density of states, as confirmed by experiment j2Hj and approximate 
theoretical techniques [24;. One thus expects that the mismatch of size produces, in 
addition to a relaxation energy E R contribution, a change in the band structure. Within 
our Augmented Space Recursion (ASR), off-diagonal disorder in the structure matrix 
because of local lattice distortions due to size mismatch of the constituents, can be 
handled on the same footing as diagonal disorder in the potential parameters [2*2] . 
The augmented space recursion with the TB-LMTO Hamiltonian coupled with orbital 
peeling allows us to compute configuration averaged pair-potentials directly, without 
resorting to any direct averaging over a finite number of configurations. In an 
earlier communication [7j we have discussed how one uses the local symmetries of the 
augmented space to reduce the Hamiltonian and carry out the recursion on a reducible 
subspace of much lower rank. If we fix the occupation of two sites, the local symmetry of 
the augmented space is lowered (this is very similar to the lowering of spherical symmetry 
to cylindrical symmetry when a preferred direction is introduced in an isotropic system). 
We may then carry out the recursion in a suitably reduced subspace. 

2.3. Static concentration wave method 

The static concentration wave (SCW) was proposed as a theory for ordering by 
Khachaturyan (|25|.|26|). The occupation probability n(f) plays the key role in this 
theory. This function n(f) that determines the distribution of solute atoms in an ordered 
phase can be represented as a superposition of concentration waves: 

n(r) = c + - E[ ( 5(^') ex P(^j " + Q*(kj) exp(-ikj ■ r)] (6) 
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where exp[ikj ■ r) is a static concentration wave, kj is a non zero wave vector defined 
in the first Brillouin zone of the disordered alloy, f is a site vector of the lattice {r}, 
index j denotes the wave vectors in the Brillouin zone, Q{kj) is static concentration 
wave amplitude and c is the atomic fraction of the alloying element. 



The study of phase stability requires accurate approximations to the configurational 
energy as well as the use of statistical models to obtain the configurational entropy. 
The configurational energy within the pair interaction can be represented in Fourier 
space as the product of the Fourier transform of the effective pair interaction V(k) and 
that of the pair correlation function Q(k): 

E^(j)^V(k)Q(k) 

where N is the number of atoms. Minimization of E will naturally occur for states of 
order characterized by maxima in the Q(k) pair correlation spectrum located in the 
regions of the absolute minima ofV(k). Consequently, much can be predicted about the 
types of ordering to be expected from a study of the shape of V{k), particularly from a 
search of its absolute minima (special points). At these points, 

\VhV(h)\ = 

This was pointed out by Lifshitz fZT\ l2*%] and Khachaturyan I2fij. Different types 
of ordered structures can be related directly to the minima of V{k). In other words, 
given the knowledge of concentration wave vectors, one can readily predict the most 
stable ordered structure of the system at low temperatures. This is comparable to the 
knowledge derived from the studies like those based on X-ray, electron and neutron 
diffraction. A peak at the T point, k = (000), indicates the phase separation, while a 
peak at the T point, k = (100), in a fee lattice suggests ordering. Peaks away from 
special points may correspond to the formation of long period superstructures. With in 
a simple mean field approximation, the instability can be obtained in the following way: 
If we add the expression for dominant quadratic term in the average energy to that of 
the configurational entropy under the simple mean field approximation we obtain an 
expression for the free energy: 

F = ^V^irii - c)(rij - c) + k B T^[rii lnn» + (1 - rn) ln(l - m)\ 

i,j i 

where Hi is the concentration of the species A at the i-th site and c is the average 
concentration of that species. If we define a configuration variable 7° as {Sn^Q (the 
symbol (• ■ -)o denotes micro-canonical averaging), which is the variable relevant to the 
stability analysis, then the harmonic term in the Taylor expansion of the above free 
energy is 

F 2 = f J2T%k(h))F(k(h))r(k(h)) (7) 

h 



Study of Phase Stability in NiPt Systems 



7 



where, k(h) = 27ih a b a and F(k(h)) = T (n(r) — c). The stability of a solid solution with 
respect to a small concentration wave of given wave vector k(h) is guaranteed as long 
as F(k(h)) is positive definite. Instability sets in when F(k(h)) vanishes i.e. 

F{k{h)) = k B T l + V(k(h)) c (1 - c) = (8) 

T l being the temperature at which the instability sets in for the considered concentration 
wave. It appears from the above expression that under a simple mean field 
approximation the spinodal is always a parabola in the (t, c) phase diagram, symmetric 
about x = 0.5. It is the concentration dependence of the effective pair interactions which 
brings about the asymmetry. 

3. Computational details 

3.1. Convergence of augmented space recursion 

The effective pair potentials are calculated at the Fermi level so one needs to be very 
careful about the convergence of the Fermi energy as well as that of the effective pair 
potentials. In fact, errors can arise in the augmented space recursion because one can 
carry out only finite number of recursion steps and then terminate the continued fraction 
using available terminators. Also one chooses a large but finite part of the augmented 
space nearest neighbour map and ignores the part of the augmented space very far from 
the starting state. This is also a source of error. 

For finding out the Fermi energy accurately, we have used the energy dependent version 
of augmented space recursion. In this version of ASR the defining Hamiltonian is recast 
into an energy dependent Hamiltonian having only diagonal disorder. We then choose 
a few seed points across the energy spectrum uniformly, carry out recursion on those 
points and spline fit the coefficients of recursion through out the whole spectrum. This 
enables us to carry out large number of recursion steps since the configuration space 
grows significantly less faster for diagonal as compared with off diagonal disorder. For 
details see ref 

We have checked the convergence of Fermi energy and effective pair potentials with 
respect to recursion steps and the number of seed energy points taking the case of NiPt3 
system. We have found that the Fermi energy and effective pair potentials converge 
beyond seven recursion steps and thirty five seed energy points. In our all calculations 
reported in the following have been carried out with eight recursion steps and thirty five 
seed energy points. 

3.2. Antiphase boundary energy 

Kanamori and Kakehasi [3U] used the method of geometrical inequalities which is capable 
of searching for ground structure. They considered the energy of the three dimensional 



Study of Phase Stability in NiPt Systems 



8 



Ising like model: 



E c — ^2 VkPk 



(9) 



k 



where, \4 is the interaction constant of the K-th nearest neighbour interaction and Pk 
is the total number of k th neighboring pairs in the given configuration. Defining the 
anti-phase boundary energy £ by 



the authors proved rigorously that for £ > Ll 2 and Ll are the corresponding 
superstructures possible at concentration 25 % and 50 % while for £ < 0, one has the 
DO22 and A 2 B 2 superstructures. We have applied these conditions in our calculations 
to find out the relative stability between DO22 and Ll 2 structures in Ni 3 Pt and NiPt3 
that between A2B2 and Llo in NiPt. 

3.3. Special-point ordering 

A wide range of phenomena related to order- disorder and magnetic transitions can be 
explained using the symmetry properties of the pair potentials (V^). If a symmetry 
element (rotation, rot at ion- inversion or mirror plane) of the space group in k-space is 
located at point h, the vector representing the gradient VhV(h) of an arbitrary potential 
energy function V(/i) at that point must lie along or within the symmetry element. If 
two or more symmetry elements intersect at point h, one must necessarily have 



since a finite magnitude vector can not lie simultaneously in intersecting straight lines 
having only a point in common. At these so-called special points, the potential energy 
function V(/i) represent an extremum regardless of the choice of the pair interaction 
energies. Thus special points play an important role in the search for lowest energy 
ordered structures. The points which differ by a vector of a reciprocal lattice are 
considered equivalent. In the case of simple structures with a single atom per unit cell, 
it is sufficient that two symmetry elements intersect at special points. These special 
points are listed in Crystallographic tables. They are always located at the surface of 
the Brillouin zone. The 'star' of a special point vector k is obtained by applying all the 
rotations and rotation-inversion of the space group on the vector k. All these vectors of 
a star are also considered equivalent. The special points of the fee structure are located 
at the points T, X, W and L of the Brillouin zone as shown in table 1. 

3.4- Ordering energy 



V 2 + AV 3 -AV A , 



(10) 



V h V(h)\=0 




The ordering energy is defined as the difference between the formation energy of ordered 
alloy and the corresponding formation energy of disordered alloy. Since we are dealing 
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k-vector Star 



Members 



Brillouin zone points 



Ordering structure 



(000) 
(100) 



1±0 



ill 

2 2 2 



[000] 

[100] [010] [001] 



1|0 §01 01| 



1|0 



111 

2 2 2 



111 

2 2 2 



§01 



OliO 





111 




2 2 2 




TTl" 




2 2 2 



r 

x 
w 



Ll 2 , Ll 
A 2 B 2 , D0 22 



Lli 



with the effective pair potentials, the ordering energy can be calculated using these pair 
potentials. The relation for ordering energy using pair potentials is given as : 

E ord= WVkSXoSXk (12) 
1 k 

where, 5x Q (Sxk) = x Q (xk) — x, x Q (xk) = 1 if the site o/k is occupied by A atom and 
x Q = if the site o/k is occupied by B atom. For Ll 2 structure (for Ni 3 Pt and NiPt 3 
in our case) the expression for ordering energy per atom in terms of pair potentials 
considering only up to fourth nearest neighbours is given as: 

Ef&Pt = -^Wi -\V2 + V,- 1 -V A ] (13) 

For LIq structure for NiPt the expression for ordering energy per atom considering up 
to fourth nearest neighbour pair potentials is given as: 

EmPt = -l\yi-V2 + V 3 -V A ] (14) 
Using these two relations we have found the ordering energy for Ni 3 Pt, NiPt and NiPt 3 . 



4. Results and discussions 



We have applied our formalism discussed in the previous section in calculating the 
effective pair potentials for the fee based NiPt alloys for concentrations x = 0.25, 0.5 
and 0.75 of Pt. The calculation of the effective pair potentials has been restricted up 
to fourth nearest neighbour interactions. Total energy density functional calculations 
were performed at the concentration x = 0.25,0.5 and 0.75 of Pt. The Kohn-Sham 
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Table 2. The calculated equilibrium lattice parameters with the choice of neutral 
charge spheres including scalar relativistic corrections. The corresponding lattice 
parameters with out relativistic corrections are given in brackets. 



concentration 
of Pt 


eq. lattice parameter 
in au. SR(NR) 


0.00(Ni) 
0.25(Ni 3 Pt) 
0.50(NiPt) 
0.75(NiPt 3 ) 

l.OO(Pt) 


6.528(6.568) (expt. 6.66) 

6.758(6.890) 

7.127(7.335) 

7.196(7.467) 
7.3685(7.683) (expt. 7.41) 



equations were solved in the local density approximation (LDA). The LDA was treated 
with in the context of linear muffin tin orbitals in the atomic sphere approximation. 
The calculations were performed non relativistically as well as scalar relativistically 
and the exchange correlation potential of Von Barth and Hedin was used. Two sets of 
calculations were performed one with the same Wigner-Seitz radius (charged spheres) 
for Ni and Pt. In other set we followed the procedure described by Kudrnovsky et al 
[TH] and used an extension of the procedure proposed by Andersen et al 20], which 
allows us flexibility in the choice of ASA radii for the constituents. The idea is to 
choose ASA radii of atomic species in such a way that the spheres are charge neutral 
on the average. The potential parameters Af and 7/ of the constituent I were then 
scaled by the factors (s 1 / s allo yy l+1 to account for the fact that the Wigner-Seitz radius 
of constituent I, s 1 , is different from that of the alloy, s alloy . These potential parameters 
were used to parameterize the alloy Hamiltonian. For the purpose of augmented space 
recursion, seven shell map was generated and thirty five seed energy point recursion 
was performed, as explained in previous section, to calculate the Fermi energy with 
the second order LMTO-ASA Hamiltonian through the recursion method using eight 
level of recursion and analytical terminator of Luchini and Nex. For the effective pair 
potentials, we used the orbital peeling method within the frame work of ASR for the 
calculation of peeled averaged Green function described in detail in the earlier section. 

In table 2, we have quoted the equilibrium lattice parameters that were used in our 
calculations. We obtained these by minimizing the total energies with respect to 
the lattice parameters. We have got slightly lower equilibrium lattice parameters as 
compared to experimental ones. This is characteristic of the local density approximation 
which overestimates bonding. 

In figure 1 we have shown the formation energy of NiPt alloy system with various 
Pt concentrations based on ordered calculations. It shows that without inclusion of 
relativistic effects the formation energy comes out to be positive which contradicts 
experimental results. With the inclusion of scalar relativistic corrections the formation 
energy comes out to be negative. This indicates that relativistic effects play an important 
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Table 3. Formation energies for Ni^Pt^with the choice of neutral charge spheres 
including scalar relativistic correction. The values in brackets are without relativistic 
correction. The corresponding estimate for charged sphere calculations are shown with 
*'s. ** refers to calculations with out combined correction. * * * refers to disordered 
formation energy 



y 


This work 


Expt.jSU 


FPLMTO+ 
fTVVMlTTl 

V V V IV J. |O^J 


LMTO 


LMTO+ 

cwmr^i 


KKR-ASA 
(KKR-CPA 


0.25 


-7.50(4.25) 
-7.59*(4.17)* 


-5.16 


-6.30 


-7.17 


-6.66 




0.50 


-9.44(4.74) 
-9.02*(4.85)* 


-7.06 


-8.69 


-8.5 


-8.95 


-12.00** pi] 

-8.io nu 

(-7.7***) PH 


0.75 


-8.15(4.22) 
-3.97*(6.65)* 


-4.78 


-6.40 


-6.70 


-9.12 





role in the stability of NiPt alloys, in agreement with earlier studies. Our results are 
in closer agreement with previous works based on the Full-Potential LMTO and the 
Connolly- Williams technique ([2l],|H2],|HHj) and with experimental estimate. Singh et 
al ^Ej have also calculated the formation energy for 50% of Pt. Their results for 
the formation energy obtained from ordered calculations without combined correction 
deviates quite a bit from ours as well as other results based on the Full-Potential LMTO 
and the Connolly- Williams technique ( [^T] , , [33J ) , which is presumably due to the 
neglect of the combined correction in reference |16j . Singh et al have also done the 
calculation with combined correction which shows better agreement. The full-potential 
methods are expected to provide better estimates than other methods. 

We next approached the problem from the disordered end. We started from a completely 
disordered alloy and set up concentration wave fluctuations in it to see when this 
destabilizes the disordered phase as suggested by Khachaturyan [23]. The calculation 
of the lattice distortion for disordered alloys has been carried out within the structural 
model given by rigid ion structure (RIS) (HSj- According to this model the lattice 
relaxes in such a way as to keep all the nearest neighbour distances close to the sum 
of the corresponding atomic radii for a particular concentration. This is found to be 
a reasonable model to deal with lattice relaxation effects in non-isochoric alloys |10j . 
Due to the distortion of the lattice, the structure matrix Sff^' (which is 9x9 matrix for 
each RR' pair and for a spd basis set) can randomly take values Sff/, Sf^/ and Sff ; 
depending upon the occupying of sites R and R'. 

Slu = S££,n R n R , + S£»[n R {l - n R .) + (1 - n R )n a \ + Sf£(l - n R ){\ - n R .) 
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Table 4. The effective pair potentials for NiPt alloy system calculated with potential 
parameters taken from calculations with the choice of charged spheres and including 
scalar relativistic correction. (O-L) refers to calculations without multipole corrections, 
M refers to calculations with multipole corrections and SCI to calculations with 
screened Coulomb interactions. US-PP refers to ultrasoft pseudo-potentials. * refers 
to non-relativistic calculations. 



Reference 


Vi 


v 2 


v 3 


v 4 




(mRy/atom) 


(mRy/atom) 


(mRy/atom) 


(mRy/atom) 


Concentration of Pt = 25% 


Present work 


11.36 


-0.05 


-0.07 


-0.41 




11.972* 


0.015* 


0.054* 


0.046* 


Concentration of Pt = 50% 


Prpepnt work 


7 832 


0.114 


-0 129 


-0 057 




8 W7* 


10* 


0^3* 


263* 


Sin cf n pi nl 1 


4.22 


1.14 


0.22 


-1 04 




4.94* 


0.52* 


0.32* 


-0.18* 


Pinski et al [H] 


9.4* 


0.8* 


0.4* 


-0.2* 


Pourovskii et al [38] 










CWM-ASA+M 


5.00 


0.25 


0.19 


-0.28 


SGPM 


5.28 


0.06 


-0.82 


-0.66 


Ruban et al [T7] 










with SGPM 










ASA+M (O-L) (SCI) 


14.05(15.44) 


0.32(-0.10) 


-1.09(-1.22) 


-1.76(-0.84) 


ASA (SCI) 


12.26(14.35) 


0.53(-0.15) 


-1.31(-1.48) 


-2.14(-0.98) 


With Connolly Williams 










ASA+M 


12.68 


1.31 


-0.02 


-0.73 


ASA+M (O-L) 


13.70 


0.49 


-0.86 


-1.39 


ASA 


14.33 


0.28 


-1.72 


-1.92 


US(PP) 


12.81 


1.30 


0.69 


-0.40 


Direct calculation(SCI) 










ASA+M 


12.45 


0.47 


-0.49 


-0.65 


Concentration of Pt = 75% 


Present work 


2.785 


0.236 


-0.116 


0.276 




3.813* 


0.361* 


-0.175* 


0.366* 
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Table 5. The effective pair potentials for NiPt alloys with potential parameters taken 
from calculations with the choice of charge neutral spheres including scalar relativistic 
corrections. The corresponding estimate for non relativistic calculations are shown 
with *'s. 



Reference 


Vi 


v 2 


V3 


v 4 




(mRy/atom) 


(mRy/atom) 


yiiisxy / aiom j 


^iiiixy/ dxoiii ) 




Concentration of Pt = 


-- 25% 




Present work 


12.34 


-0.092 


-0.046 


-0.54 




13.08* 


-0.021* 


U. ioz* 


-U.U41* 




Concentration of Pt = 


= 51) /o 




Present work 


10.08 


0.1 


c\ c\c\ a 

0.004 


-0.24 




10.111* 


0.126* 


U.Z4D* 


(J.Li 0* 


Singh et al [To] 


16.02 


1.34 


0.06 


-1.58 




11.96* 


0.66* 


0.28* 


-0.46* 


Ruban et al [TT] 










Neutral(GPM) 


5.49 


1.22 


0.01 


-0.73 




Concentration of Pt = 


= 75% 




Present work 


8.9 


0.26 


0.1 


0.02 




7.874* 


0.297* 


0.276* 


0.34* 



where 




1 if R is occupied by A 
if R is occupied by B 



Considering the example of calculation of Sf£ where B is the larger atom (e.g. Pt in 
the present case), this matrix for a specific pair among 12 nearest neighbours connects 
an A atom at the site (0, 0, 0) and a B atom, which in the undistorted case would have 
been at the position (| , | , 0) is now at ((| + d), (| + d), d), where d is the displacement 
due to lattice distortion and a is the lattice constant. We have assumed that the lattice 
expands equally in the x, y and z directions. With these new coordinates and assuming 
that all other neighbouring coordinates are fixed at undistorted fee positions (which 
is the essence of terminal point approximation JH] ), we have computed the structure 
matrices Sf^ h SfE/ an d SEE- This takes into account both the effect of radial distortion 
as well as angular distortion (the nearest neighbour is now \J ^ + 2a<i + 3d 2 instead of 
-7= and the nearest neighbour vector is ((f + d) , (| + d) ,d) instead of (f , | , 0) in the 
above example). The values of d for Sf£ came out to be 0.064 a, 0.052 a and 0.054 a 
for 25%, 50% and 75% concentration of Pt. The details of the calculation scheme can 
be found in reference [THj. In figure 2 we have shown the relative magnitudes of nearest 
neighbour distances for different concentrations of Pt in NiPt alloy system compared 
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Table 6. The anti-phase boundary energies for Ni^Pt^ alloys from charged and neutral 
sphere calculations 



concentration 
of Pt 


APB energy 


(mRy/atom) 


Charged spheres 


Neutral spheres 


Relativistic(Non-Relativistic) 


0.25 


1.41(0.017) 


2.07(0.793) 


0.50 


-0.402(-0.94) 


0.876(0.158) 


0.75 


-1.804(-2.525) 


0.06(-0.553) 



to Vegard's law values for average bond length. We have computed the effective pair 
potentials for two sets of potential parameters with charged and charge neutral spheres. 
Figure 3 shows that the effective pair potentials for NiPt 3 is very small in magnitude 
using potential parameters with charged spheres. We even used these pair potentials 
and calculated the anti-phase boundary energy according to the prescription described 
in the previous section. The anti-phase boundary energy comes out to be negative for 
NiPt3 and NiPt indicating stability of DO22 over Ll 2 for NiPt3 and A 2 B 2 over Ll 
for NiPt. Further we calculated the minima of the special points according to the 
prescription described in the previous section. In the case of NiPt 3 and NiPt shown 
in figure 5, we could not get the minima at (100) which is not quite correct because 
experiments show NiPt 3 has Ll 2 and NiPt has Ll ordering. But in the case of Ni 3 Pt we 
could get the positive anti-phase boundary energy as well as minima at (100) correctly 
showing the ordering Ll 2 . In figure 4 we have plotted the effective pair potentials as 
a function of energy relative to Fermi energy and number of neighboring shells with 
charge neutral potential parameters including scalar relativistic correction which shows 
that the first nearest neighbor pair potentials are larger in magnitude than the second, 
third and fourth nearest neighbour pair potentials. With potential parameters from 
neutral sphere calculations including scalar relativistic correction for NiPt 3 and NiPt 
the anti-phase boundary energies come out to be positive and the minima of special 
points are at (100) correctly showing Ll 2 and Ll orderings. If we use charge neutral 
potential parameters with out including scalar relativistic effect the anti-phase boundary 
energies come out to be positive for Ni 3 Pt and NiPt but negative for NiPt 3 . This shows 
for NiPt 3 both scalar relativistic as well as charge transfer effects play important role 
to predict correct ground state. 

So, we argue that on increasing the concentration of Pt atom the careful treatment to 
take into account of charge transfer effect becomes increasingly important. In figure 4, 
we have also shown the effective pair potentials without scalar relativistic corrections. 
For NiPt 3 it is clearly seen that the effective pair potentials with scalar relativistic 
correction are larger in magnitude than the non relativistic ones which is expected 
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because of higher concentration of Pt. 

In figure 6 we have plotted that the effective pair potentials vs concentration of Pt 
with charge neutral potential parameters including scalar relativistic correction which 
shows that the first nearest neighbour effective pair potentials decrease with the increase 
of the Pt concentration. Singh et al ^Hl EH] have also calculated the effective pair 
potentials using KKR-CPA-GPM method. Their values of effective pair potentials are 
much larger than ours. They pointed out that due to the large values of effective pair 
potentials the ordering energy and ordering temperatures (transition temperatures) are 
much higher than that observed experimentally. Our estimates give rise to instability 
temperatures which are closer to the experimental results (shown in figure 7). For 
example, our estimate for the instability temperature for the 50% alloy is 1683°K, 
whereas the estimate from the KKR-CPA is around 2979°K. The experimental estimates 
of the transition temperature is 950°K |SZj. In KKR-CPA-GPM method one considers 
only the single site approximation and one does not take into account any off diagonal 
disorder which may arise because of size mismatch of the constituent atoms. The ASR, 
on the other hand, as discussed earlier can do this with facility. Our test calculation for 
NiPt (50% concentration of Pt) without taking into account lattice relaxation due to size 
mismatch effect gives an estimate of instability temperature of 2363°K which is indeed 
higher than that of our original estimate with taking into account lattice relaxation due 
to size mismatch effect. Furthermore, Singh et al pHlEE] in their calculation for charge 
neutrality have taken the ratio of Wigner Seitz radii of Ni and Pt as 0.95. We, on the 
other hand, have varied the ratio, with the provision that the total volume is conserved, 
till, on the average, the spheres become charge neutral. We have observed the ratio 
to be 0.909, 0.913 and 0.919 for the Ni 3 Pt, NiPt and NiPt 3 . Given these calculational 
differences, it is not surprising that our calculations result in smaller values of the pair- 
potentials leading to better estimates of the instability temperatures. The calculations 
of Pinski et al J21 were carried out without scalar relativistic effects. Their values are 
consequently rather large as compared to ours. Ruban et al [T^j have calculated pair 
potentials for 50% concentration of Pt using different methods and showed that different 
methods give different values of pair potentials. Their nearest neighbor pair-potential is 
slightly higher than ours. The effective pair potentials obtained by Pourovskii et al [3B] 
from the neutral charge spheres GPM method are similar to the estimates of Ruban et 
al 

In figure 7 we have shown the ordering energy, anti-phase boundary energy and 
instability temperatures vs concentration of Pt with charge neutral potential parameters 
including scalar relativistic correction. The ordering energy in all three cases NisPt, 
NiPt and NiPt 3 is negative showing the stability of ordered structures compared to 
disordered solution. Among all three concentrations, NiPt attains maximum value of 
ordering energy which confirms that Ll in NiPt system is the most stable structure. 
The anti-phase boundary energy in all these cases NisPt, NiPt and NiPt3 comes out 
to be positive showing the ordering structures LI2 for NisPt, L1q for NiPt and LI2 
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for NiPt.3 as described above. The magnitude of instability temperatures using the 
charge neutral potential parameters comes out to be larger than the experimental 
transition temperatures. However, the qualitative trend of the change of the instability 
temperatures with changing concentration of Pt is right. The calculated qualitative 
phase diagram (instability temperature vs concentration of Pt) shows asymmetric 
feature which is not observed experimentally. This could be due to the neglect of 
magnetism in the calculations of effective pair interactions which can have significant 
effect particularly in the high concentrations of Ni. Amador et al jHI] also reported 
the phase diagram (instability temperature vs concentration of Pt) of this system 
described by the nearest neighbour tetrahedron effective interactions from clusters with 
appropriate effective volume. Their values for transition temperatures are smaller than 
ours but there is high asymmetry in their phase diagram and even the trend is not same 
as experimental findings. 

5. Conclusion 

Our total energy calculations for the ordered alloys indicate that in order to have the 
correct sign for the formation energy, it is essential to include relativistic corrections. 
Our analysis of the concentration wave approach indicates that for Ni 3 Pt neither 
relativistic correction nor the charge transfer effect is essential for the correct prediction 
of the LI2 ground state. For NiPt although scalar relativistic correction is not essential, 
careful treatment of charge transfer effect is a must to predict the correct ground state 
(LIq). For NiPt3 both these corrections are essential to predict the correct ground state 
Ll 2 . 

Although it seems that qualitatively the relativistic corrections and charge transfer effect 
plays the essential role only for the high Pt content alloys, for quantitative prediction of 
the instability temperature both these corrections are required across the concentration 
range. 

The main conclusions of this paper are : 

• We have demonstrated that for accurate prediction of the ground state structures 
and instability temperatures for alloys with components with large atomic size 
differences like NiPt, it is essential to take into account both relativistic corrections 
and averaged charge neutrality of the atomic spheres. 

• We have also demonstrated the augmented space recursion combined with the first- 
principles tight-binding linearized muffin-tin orbitals and the orbital peeling are 
both computationally feasible and suitable techniques for such studies as described 
above. 

These techniques will form the basis of our further study into similar alloy systems, but 
with magnetic effects included. 
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Figure 1. Formation energy vs concentration of Pt with the choice of neutral charge 
spheres 
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Concentration of Pt 



Figure 2. Nearest neighbour distance vs concentration of Pt with the choice of neutral 
charge spheres. For comparison the average bond length given by Vegard's law is shown 
in dashed line. 
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Figure 3. (i) The effective pair potentials with potential parameters taken from 
calculations with the choice of charged spheres including scalar relativistic corrections. 
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Figure 4. (i) The effective pair potentials as a function of energy relative to 
Fermi energy with charge neutral potential parameters including scalar relativistic 
corrections, (ii) Comparison between the first nearest neighbour effective pair potentials 
with scalar relativistic corrections and without scalar relativistic corrections by taking 
charge neutral potential parameters, (iii) The effective pair potentials as a function 
of shell numbers with charge neutral potential parameters including scalar relativistic 
correction. 




Figure 5. The V (k) surface for NiPt alloy system with potential parameters 
calculated with the choice of (i) charged spheres (ii) charge neutral spheres, on k z 
= plane. The figures in inset are corresponding rescaled V (k) surfaces on k z = 
plane along the (100) to (110) direction to view the minima. 
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Figure 6. The effective pair potentials vs concentration of Pt with the choice of 
potential parameters with (i) charged spheres including scalar relativistic correction 
and (ii) charge neutral spheres including scalar relativistic correction. 
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Figure 7. Ordering energy, anti-phase boundary energy and instability temperatures 
vs concentration of Pt with the choice of charge neutral potential parameters including 
scalar relativistic correction. 



